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Abstract 

Gross-Pitaevskii equation for Bose-Einstein condensate confined in elongated cigar-shaped trap 
is reduced to an effective system of nonlinear equations depending on only one space coordinate 
along the trap axis. The radial distribution of the condensate density and its radial velocity 
are approximated by Gaussian functions with real and imaginary exponents, respectively, with 
parameters depending on the axial coordinate and time. The effective one-dimensional system is 
applied to description of the ground state of the condensate, to dark and bright solitons, to the 
sound and radial compression waves propagating in a dense condensate, and to weakly nonlinear 
waves in repulsive condensate. In the low density limit our results reproduce the known formulae. 
In the high density case our description of solitons goes beyond the standard approach based on 
the nonlinear Schrodinger equation. The dispersion relations for the sound and radial compression 
waves are obtained in a wide region of values of the condensate density. Korteweg-de Vries equation 
for weakly nonlinear waves is derived and existence of bright solitons on a constant background is 
predicted for the dense enough condensate with repulsive interaction between the atoms. 
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I. INTRODUCTION 

Formation of Bose-Einstein condensate (BEC) in dilute gases 

HQS 

has created a 

new field of nonlinear physics with new specific problems. Among new features of BEC in 
gases are the spatial inhomogeneity created by the confining potential, interplay of coherent 
and nonlinear phenomena, multi-species condensation, etc. Many of such new problems 
have already been addressed in current literature. In particular, the condensate oscillations 
near the center of the trap and its free expansion after switching off the trap potential was 
considered in Refs. ll| . Creation of strongly elongated ( "cigar-shaped" ) 



traps has made possible to generate dark 
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13| and bright 3, |15| sohtons in BEC. 



However, theoretical description of the solitons proves to be not a simple problem. Indeed, 
the standard reduction of the Gross-Pitaevskii (CP) equation to the one-dimensional (ID) 
nonlinear Schrodinger (NLS) equation is justified only under the condition that axial size of 
the solution (along the trap axis) is much greater than its radial size. However, the actual 
solitons observed in experiment Q, Q| have comparable axial and radial sizes. An attempt 
to go beyond the NLS approximation was made in Ref. under the assumption that 
variation of the field variables in the axial direction takes place on much longer scale than 
that in the radial direction. Therefore the equation with a "non-polynomial interaction" 
obtained in Ref. [igI describes the ground state of BEC with repulsive interaction when the 
above assumption is correct but its application to "short" dark or bright solitons is not well 
justified. 

In this paper we consider dynamics of BEC in the cigar-shaped trap by means of vari- 
ational approach without assumption that the field variables change slowly along the trap 
axis. The radial distribution of BEC density and its radial velocity are approximated by real 
and imaginary Gaussian functions, respectively, with parameters depending on the axial co- 
ordinate z. This approach allows us to reduce the three-dimensional (3D) BEC Lagrangian 
to an effective ID one which depends on the axial wave function \I^(z,t), the mean con- 
densate radius b{z,t), and the radial velocity potential a{z,t). The resulting system of 
equations for the above three variables, though being rather complicated, has an advantage 
of dependence on only one spatial coordinate z. It contains as limiting cases the NLS and 
the "non-polynomial interaction" approximations. The derived effective system allows one 
to investigate bright and dark solitons in a wide domain of their existence. In particular. 



2 



the bright sohtons can be considered right up to the collapse instability threshold. We also 
study the linear waves propagating along the condensate with repulsive interaction between 
the atoms. There are two types of waves: the low frequency sound waves and the high fre- 
quency radial compression wave. The dispersion relations for the linear waves, applicable to 
the high density BEC, are obtained in the long wavelength limit. In the next approximation 
of "small but finite" amplitude of the ground state perturbation we obtain the Korteweg-de 
Vries (KdV) equation for weakly dispersive and nonlinear waves propagating along the con- 
densate. An interesting feature of this equation is that in a high density BEC the coefficient 
in the dispersion term changes sign. An immediate consequence is that fast solitons of small 
amplitude, propagating in a repulsive BEC, are bright rather than dark, contrary to what 
is usually assumed by analogy between the CP and ID NLS equations. 



II. GENERAL QUASI-ID EQUATIONS FOR BEC EVOLUTION 

The starting point for our consideration is the CP equation for the wave function i/j of 
the condensate, 

n'^ 1 

ihipt = A?/" + -mulr'^ij + V(z)ip + glipl^ij, (1) 

2m 2 



where r = a/o;^ + is the radial coordinate, g = Aixh^as/m is the effective nonlinear 
coupling constant, denotes the s-wave scattering length, m is the atomic mass, uo± is the 
frequency of atomic oscillations in the radial direction, and V{z) is the trap potential in the 
axial direction of a cigar-shaped trap. The condensate wave function is normalized to the 
number N of atoms, 

j iV'prfr = A^. (2) 

The CP equation admits the least action principle formulation with the action functional 

S = j Ldt, L = J C- 2nrdrdz, (3) 

written under assumption that the condensate motion is axially symmetric, where the La- 
grangian density is given by 

>C = y m - r^t) + l^lVV^r + ^mcuVl^r + VizMl' + (4) 

It is well known (see, e.g. that if the linear (ID) density rii = \ip\'^27irdr of the 
condensate is small enough, 

\as\ni < 1, (5) 



then the 3D GP equation can be reduced to the ID NLS equation 



2m 



where 



giD 
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,2 ' 



n 



(6) 



(7) 



27ra]_ ma\ V "itluj^ 

This reduction can be performed, for example, by variational method, when under suppo- 
sition (jSJ the transversal degrees of freedom are "frozen" and transverse motion is reduced 
to the ground state oscillations in the radial direction. As a result, the condensate wave 
function ■?/' can be factorized as 



exp 



2a\ 



(8) 



Its substitution into 0, © and integration over radial coordinate yield an effective ID 
Lagrangian and the corresponding Lagrange equation for '^{z,^) is Eq. (jH} (see, e.g. 
Here we want to generalize this procedure to the case when la^jni ~ 1 and the radial motion 
of the condensate must be taken into account. 

Assuming quasi-lD geometry, we approximate the wave function as 



\/7Tb{z,t) 



exp 



exp ( -a{z,ty ) 



(9) 



2b^{z,t) 

Here b{z, t) characterizes the local radius of the condensate and a{z^ t)r gives the local radial 
velocity. Substitution of into (jH) and the result into © with integration over the radial 
coordinate yield an effective ID Lagrangian 



ih 



n 



2m 



1*1 



h 



^ + aV + ^ + ^a'X 1 1*1" + -a,h' (**! - 



,21,4 



2m \b^ b^ 2 

To simplify the notation, let us introduce dimensionless variables: 

u 



(10) 



r = Lj±_t, z = a±(, 



a_|_w, a 



a I a± 



(11) 



In the dimensionless variables the ID Lagrangian becomes 



ID 



ma , 



1 



u\ 



(12) 
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Thus, in the above reduction, the evolution of BEC is described by the complex longitudinal 
wave function u{(,t) and two real functions, w{(,t) and t), corresponding to the local 
mean radius of the condensate and its radial velocity, respectively. Equations, governing the 
evolution of these variables are to be obtained from the action principle 



S 



Ldt = min, 



with the effective Lagrangian density given by Eq. (jl2p . We have: 



1 



iur + -ua - y{C)u - 3(7 „ 
2 



W 



W 



w 



i\u\ 



)c ^2„.,4 ^ a „..2 



\U\ 



(13) 



(14) 



(15) 



|up ^ 2|Mp 

where the longitudinal wave function is normalized as follows 



\u\^dC = N. (17) 

Thus, we have transformed the 3D GP equation to a quasi- ID form for the case of elongated 
cigar-shaped geometry when the radial distributions of the condensate density and its radial 
velocity can be approximated by simple Gaussian functions. The variables in Eqs. p4|l - 
(I16|) depend only on one spatial coordinate — a strong advantage in numerical simulations. 
Besides that, system fjl4|) - (fTm) can be analyzed analytically in some important limiting cases 
(see below). 

For some applications it is convenient to cast Eqs. ()14|) - |T6|) into a hydrodynamic-like 
form by means of the substitution 



u 



^exp(^z I v{C,T)dC^ , (18) 



where p{(,t) is the density of BEC and v{(,t) denotes its axial velocity. As a result, we 
obtain the hydrodynamic system: 

p.+ (p(t; + X/3c)), = 0, (19) 



P j , 1/ Pec , 2 , w^C , ^CPC 1 ,p2„„4 
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rT-^ + ^ + — + ^-o/^^^ =0, (20) 



p,^^.^ + "£i£L + L l.0^.f^^ + ,^^ + 2G-L, (21) 



{w^p)r + {w^pv + wV/3c)c = 2/5wV. (22) 

In the following sections we apply the general quasi-lD equations to some important par- 
ticular cases. 



III. STATIONARY SOLUTIONS 

In a stationary state all velocities are equal to zero (t; = 0, /5 = 0), density p does 
not depend on r and the condensate wave function m(C, t) depends on r only through the 
phase factor exp(— ipr), where p is dimensionless chemical potential. Setting /3 = and 
introducing the stationary variables f/(C) and o"(C) by the formulae 

«(C, r) = e-^^-f/(C), w\C, r) = l/a(C), (23) 
after simple algebra we obtain from Eqs. (|n|l - (|T6|) the following system 

% + - 2^(0 - - = 0, (24) 



2 

- 2GaU^ = 0. (25) 



2 \ a J ^ \ a 
System (f^ - ipHjl has the energy functional 

£ = 



1 erf 



dC- (26) 



Here integrand can be considered as a Lagrangian with the spatial variable ( playing the 
role of time. It is worth noticing that if we introduce instead of a an angle 6 by the relation 

cT = exp(2^), (27) 
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then the energy functional takes the form 



- {Ul + U^el) + (y{C) + cosh 20 - ii))U^ + Ge^'U^ 



dC 



(28) 



_2 v-'C 

which is formally an action functional of a particle moving on the plane, where the couple 
{U, 6) assumes the role of the polar coordinates and the external potential is 

$(f/, e) = {ii- V{C) - cosh 2e)U^ - Ge^^U\ (29) 



A. Small-density limit 

First of all let us consider the limit of a low density BEC: 

|G|f/2 < 1. (30) 

In the dimensional variables, with account of the estimate [/^ ~ |^pct± ~ nia± (where 
rii is the longitudinal density of the condensate), condition (fHUjl coincides with Since 
the characteristic length in the axial direction cannot be less than unity (i.e. a± in the 
dimensional units), we find in this limit from ()25p 

a = l + (Ti, ~ |G|f/2 < 1. (31) 

Taking estimate (|^ into account, we reduce Eq. to the standard (stationary) ID NLS 
equation 

% + 2(/i - 1 - V{C))U - 4Gf/^ = 0. (32) 
This equation is studied very well and we shall not consider it any further. 



B. Ground state of BEC with repulsive interaction {G > 0) 

Another important limiting case corresponds to very slow dependence of BEC parameters 
on the axial coordinate C, so that one can neglect the derivatives in Eqs. (j21I), ()25|) . If we 
denote by Z the characteristic axial length, then the condition for neglecting the spatial 
derivatives can be written as 

Z"^ < G'af/^ (33) 
or in the dimensional units, with Z ~ //a_L, G = as/aj_, U"^ ~ n\a_[_ ~ Na±/l, cr ~ 1, as 



If condition ()34|) is satisfied, then Eq. ()25|) gives the relation 

a = (1 + 2GU^)-^'^ (35) 

equivalent to one obtained in Ref. [l6]. Substitution of into and dropping off the 
derivatives relates the field U to the external trap: 

Solution of this equation with respect to U'^ yields the condensate density as a function of 
the axial coordinate 

P(C) = = ^ ^ - ^(0 + V(/^-nC))^ + 3]'- l} . (37) 

This formula corresponds to well-known Thomas- Fermi approximation for condensate's den- 
sity distribution. The density p vanishes at C = where value of is determined by the 
equation 

/i - = 1 (38) 

which gives the relationship between p and Zm- The complementary relation between Z^ 
and IX is given by the condition 

j p{C)dC = N, (39) 

where the integration is taken over region where p > 0. These two relations allow us to 
express /i and Z^ as functions of the number of atoms N. 

In the low density limit, i.e. GU^ ^ 1, Eq. reduces to the usual formula 

P=^(/^'-nC)), /i' = p-l, (40) 

which also follows from the ID NLS equation ()32j] of the preceding subsection. 
For the parabolic potential 

V{0 = 1^\', A = ^, (41) 

A UJ_\_ 

from Eq. ^ we obtain p' = ^A^Cm so that Eq. ^ with account of ^ yields the well- 
known expression for density distribution in the Thomas-Fermi approximation: 

PiC) = U' = ^(l-^], Cm = {3GN\-Y^\ (42) 



AG V 
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C. Dark solitons in BEC with repulsive interaction 



It is well known that the NLS equation ()32j) with V{() = and G > has a dark soliton 
solution 

f/(C) = Uo tanh{UoV2G C) (43) 

under condition ()30|) . that is fi' ~ GUq <^ 1. The soliton width is much greater than unity, 
1/{Uo^/G) ^ 1, and it must be much less than the longitudinal size (m of the condensate 
estimated in the preceding subsection. For BEC confined by a parabolic trap we have 
Uq ~ N/ (rn SO that the formulated condition reads 

1 GN , ,,,, 
« « 1, (44) 

Sm S-m 

where the left inequality actually coincides with condition ()33|) in the small density limit 
and the right inequality is equivalent to (jHUI) . Substitution of Cm from Eq. yields the 
applicability condition of the soliton solution ()43|) in the form 



y/X<^GN <^ X-\ (45) 

For large enough the right inequality in Eq. is violated and we have to return 

to the system (j^^ and ()25|) . Since the left inequality in Eq. is easily fulfilled for a 
dense BEC in a cigar-shaped trap, we suppose that the longitudinal size ~ Cm of the whole 
condensate is much greater than the soliton width. To find the soliton solution we take 
V{() = and require that U satisfies the following asymptotic boundary conditions: 

\U\^Uo at C ^ ±oo. (46) 
Then Eq. (17^ (see also Eq. (jHH)) gives 

a^ao = (l + 2G'f/o^)-^/2 ( ^ ±oo, (47) 

while from Eq. pljl (see Eq. (j36|l ) we find the value of /i: 

u = ^ + ^^^0 (48) 

Thus we arrive at the system: 

% + [2/i - a - i - ^] t/ - AGaU' = 0, (49) 




1 (a^U^ 
2 



a 



C 



a - - ] - 2GaU^ = 



whose solution minimizes the energy 

'1 



£ 



-Ul + GaU' + l(a + -+ 



(50) 



(51) 



2y ' a 4a2^ 

where /z is fixed by the asymptotic value Uq of the wave function as dictated by formula 

iHl). 

The solution of system was found numerically and the results are presented 

in Fig. 1. As we see, even for relatively small values of Uq (scaled by \/G in the figure) 
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FIG. 2: The number of atoms in the dark sohton "hole" as function of the asymptotic vahie GUq. 
The solid line corresponds to our quasi-lD approach and the dashed line to the NLS approximation. 

the transverse radius of the density distribution changes quite considerably along the trap 
axis rather than being constant as is supposed in the NLS equation approximation. The 
soliton width is about unity (in the dimensionless units) for our numerical solutions which 
also contradicts the condition of applicability of the NLS equation. 

The discrepancy between the improved quasi-lD description and the NLS solution can 
be manifested by comparison of the number of BEC atoms expelled from the dark soliton 
"hole", defined as 



In Fig. 2 we give the plot of the number of expelled atoms vs. GUq. It shows very convinc- 
ingly that the NLS approximation is quantitatively incorrect outside the region GUq ^ 1. 

We conclude this section by noticing that our quasi- ID theory permits one to describe 
dark solitons in the high density BEC when the standard NLS approximation breaks down. 

D. Bright solitons in BEC with attractive interaction (G < 0) 

In the small density limit (j3(Jj) the condensate with attractive interaction between the 
atoms is also described by the NLS equation (jH^ . now with G < 0. In this case BEC can 
be confined in the axial direction of the cigar-shaped trap by the interatomic interaction 




(52) 
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only. Supposing the latter, we can omit the longitudinal trap potential V{() and obtain the 
bright soliton solution 

[/(C) = , (53) 

cosh([/ov^C) 

where Uq is the amplitude of the soliton connected with the number of atoms N = J U'^d( 
by the relation 

(54) 




and /i' = fi — 1 = —\G\Uq. Then condition (pUj) of applicabihty of the NLS equation can be 
written in the form 

l^'l = \G\U^ ~ i\G\Nf < 1. (55) 

It is clear that this condition breaks down for large enough number of atoms ~ Vl^l 
one cannot neglect the effect of the atomic interaction on the transverse size of BEC. The 



transverse degrees of freedom lead to the collapse instability of BEC (see, e.g. IL^), which, 
as is known, also takes place in BEC confined in a cigar-shaped trap. Numerical results of 
Ref. jiof indicate that BEC in the cigar-shaped trap collapses for 

|G|A^ > 0.676 (56) 

and in this region of parameters the soliton solution ceases to exist. Our quasi-lD approach 
captures this essential property of the attractive BEC, and hence, in our approach, the 
bright soliton solution can be studied in the whole region of its existence. 

We have solved numerically the system (j49j) - (j5(jp with G < under the boundary condi- 
tions 

U{C) ^0 at C ^ ±oo. (57) 

In Fig. 3 we show the wave function amplitude and transverse radius profiles of BEC for 
several values of \G\N. It is clearly seen that the amplitude of the soliton increases while 
its transverse radius decreases with growth of IGjA^. For small \G\N the solution is well 
approximated by the NLS solution (jSSI), especially the amplitude profile. At C ^ ±cxd the 
radius approaches 1, what corresponds to the small amplitude limit when the transverse 
wave function is given by the oscillator ground state. Fast decrease of w in the center of the 
soliton with growth of \G\N precedes the collapse instability of BEC. 
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10 15 



FIG. 3: The wave function amplitude (a) and transverse radius (b) of BEG corresponding to the 
bright sohton solution vs. C for several values of |G|A^. 



To clarify the transition to collapse, in Fig. 4 we show the dependence of /i on IGjA^ in 
our quasi-lD approximation (the solid line) and in the NLS equation case (the dashed line). 
We see that dN/d^ vanishes at 

|G|A^ = 0.73 (58) 

which is in reasonably good agreement with the exact 3D critical value given by 

Finally, the axial width of our soliton solution is about unity (i.e. ~ a± in the dimensional 
units) in qualitative agreement with the features of bright solitons observed experimentally 



in Refs. 



1$. 
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FIG. 4: The dependence of /i on |G|A^ in quasi-lD approximation (solid line) and in the NLS 
equation case (dashed line). 

IV. LINEAR WAVES IN CONDENSATE WITH REPULSIVE INTERACTION 

Now let us consider the linear waves propagating along a cylindrical condensate, that is 
we assume that the wavelength is much less than the axial size of the condensate so that 
the influence of the axial trap can be neglected. 

We shall start from equations in the hydrodynamical form (fTI]|) - (f^^ with V'{() = in 
Eq. ()20j). In the stationary state the condensate has constant linear density po a^nd the 
transverse radius wq related with po by the equation (see (123)) 



In the presence of a linear wave p and w are slightly deviated from these stationary values, 
and a wave can be described by the following small variables 



Linearization of the system p9|) - ()22|) with respect to these small variables leads to the 
following equations: 



= 1 + 2Gpo. 



(59) 



p' = p - Po, w' = W - Wo, V, p. 



(60) 



(61) 




(62) 
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w'^ - wo(3 + iu^o^/^cc = 0, (64) 

We look for the solution in the form (p', w', v, (3) oc exp[i(irC — l^r)], the most interesting 
limiting cases of which can easily be investigated. Indeed, in the long wavelength limit K 
we can neglect the terms with higher derivatives. Then Eq. (jMj) gives (3 = w'^/wq and the 
substitution of this into Eq. yields the equation 



w. 



+ 4w;' = 0. (65) 



TT 3 

For low frequency sound wave ^ 4 and we obtain w' = {G/2wq)p' . Then after elimination 
of the radial variables (3 and w' the system (j6H) - (|63|) reduces to a couple of equations for 
the longitudinal variables: 

pt.Po.c^O, . 1 . . 0. (66) 

These lead at once to the expression for the sound velocity 

2 ^ Gpo{2 + 3Gpo) 
K' (l + 2Gpo)'/'' ^ ' 

In dimensional units we have 

Gpo = ciswlni (68) 

where 

ni = 'na\n (69) 

is the one-dimensional density in the limit of tight radial trapping (Gpo ^ 1) and n is the 
three-dimensional density evaluated in the "center of the trap" at r = 0. Thus Wq is obtained 
from the equation Wq = 1 + 2awlni which yields 



Wo 



1/2 



a/1 + (ani)2 + ani . (70) 



Then, in the dimensional units, we obtain the following expression for the velocity of sound 
propagating along the trap: 

.2 / ^ V 2 + 3Gpo 
15 



In the low density (tight radial trapping) limit we have Gpn ~ asUi -C 1, so that formula 
(f7T|) reproduces the well-known result (see, e.g. Refs. [3,01) 

h 



-y/2ani. 



(72) 



Expression (f7T|) gives dependence of the sound velocity on the density ni correct up to 
Gpo ~ 1. 

Now, if we neglect all (^-derivatives in Eqs. ()6H) -()63 p . then we obtain p' = 0, v = 0, i.e. 
there is no motion of the gas along the trap, and Eq. leads to the frequency 



n = 2 



(73) 



of compressional oscillations in radial direction. This formula agrees with the general state- 



ment ut 



21 



|22| that radial oscillations of a very elongated condensate has the frequency 
Q = 2u!± (in dimensional units) independent of interatomic interaction. 

Equations (jUT|) - (jn^ permit us to find corrections to the frequencies in equations (|U7jl and 
(I73p . We have the following dispersion relation for the sound waves 



« 4 



(Gpo) 



3(G'po 



2Gpo)= 



It is worth noticing that the dispersion correction changes its sign at GpQ 
order dispersion correction has to be taken into account at this point. 
For high-frequency radial compression wave we obtain 

2 + 6Gpo + 5(Gpo)' 



(74) 

3.91 and higher 



4 + 



(l + 2Gpo)3/2 

Naturally, a correction term is not universal and depends on the atomic interaction. 



(75) 



V. WEAKLY DISPERSIVE AND NONLINEAR WAVES IN REPULSIVE EEC 

At last, let us turn to waves of small but finite amplitude which propagate in repulsive 
BEC confined in an elongated cigar-shaped trap. 

As we know from the preceding section, the linear waves of long wavelength propagate 
with the sound velocity Cg. For finite wavelength we have to take into account the dispersion 
correction described by the second term in the dispersion relation (f7^ . We must also take 
into account small correction arising due to nonlinear terms in the system (|T ^ -(|^ ^ (with 
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V'{() = 0). To this end, we use the standard singular perturbation theory (see, e.g., 23|) 
and introduce the stretched variables 

T = e^/\ Z = e'/\C + CsT) (76) 

corresponding to the slow evolution of the wave in the reference system moving with the 
sound velocity Cg. Here e is a small parameter which characterizes the magnitude of a 
perturbation of the ground state due to presence of wave. Hence, we introduce the series 
expansions in powers of e: 

p = po + ep^^^ + eV^^^ + • • • , w = wq + ew^^^ + e^w^^^ + . . . , 

y = e^(i) + + . . . , (3 = ^3/2^(1) + p^/^f^m + . . . . ^^^^ 

Substitution of these series expansions into Eqs. (fT ^ -(P ^ gives, in the first two orders in 
powers of e the ground state formula (j^Uj) . the relationships 

^« = ^P«, .« = -V^), /?« = ^P«, (78) 
2w^ Po 2w^ 

and reproduces the formula (j67|) for the sound velocity. The last relevant terms of the 
expansions in powers of e yield the equations 

Wq Wq Wq Wq Wq 

P? + + ipot"?4'l = -<:,pf - Po4". (80) 



_ _ (2) 2(3G'po + 1) (2) _ 3G (2) 

- CgV^ + 3 Wz oPz ' 

Wq Wq 

which can be considered as a linear system for the variables p^^\ w^^\ v""^ (we have omitted 
the equation following from Eq. (j^^ which gives the expression for in terms of the 
other variables). It is easy to find that the right-hand sides of Eqs. ()79 |) -()81 |) are linearly 
dependent. Hence, for existence of a non-trivial solution, the left-hand sides must satisfy the 
corresponding compatibility condition which with account of Eqs. (|7F?|l can be transformed to 
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the KdV equation taking after transition to the "laboratory" reference frame with variables 
r, ( and the density disturbance p' = ep'-^^ the form: 



8c. 



^ , (Gpo)^(l-3(Gpo)^) ] , 3G , (Gpo) 



;i + 2Gpo) 







2(l + 2Gpo) 



pVc = 0. (82) 



The linear terms here reproduce the first two terms (ff^ of the expansion of the dispersion 

relation in powers of K. This KdV equation has well-known soliton solutions. 

In the low density limit Gpo ^ 1 Eq. (|82j) takes the form 

1 3G 

pt - csp'^ + —p'^^^ p'p'^ = 0, Gpo « 1. (83) 

It is worth while to note that this equation can be derived directly from the ID NLS equation 
(jni) which is a correct approximation to the GP equation in this limit. Indeed, evolution 
of a small-amplitude disturbance a(z, t) of the ground state \I^o = const, i.e. \I' = [\I'o -l- 
a{z^ t)] exp(— itut + ^0(2:, t)) is governed by the KdV equation (see, e.g. [2J|) 

at - Csa^ + a^zz r — aa^ = 0, (84) 

where Cg is the dimensional sound velocity ()72|) . Then transformation of this equation 
to dimensionless units with account of the relationship p' = 2a±'^Qa between small 
disturbances of the density (p') and the amplitude (a) reproduces Eq. (|H^. 

A remarkable feature of Eq. (jH^ is the fact that the dispersion term changes sign at 
Gpo = 3.91 while the coefficient at the nonlinear term is always negative. This means that 
the small amplitude solitons propagating on a constant background in repulsive BEG are 
dark only in the region GpQ < 3.91 but become bright in the region Gpo > 3.91. Obviously, 
such behavior is explained physically by competition of two mechanisms of energy changes 
due to presence of the wave — one is caused by the change of the potential energy of the gas 
in the transverse trap and another one is related with the nonlinear energy of interatomic 
gas interaction. Thus, we have found that bright solitons can propagate along dense enough 
repulsive BEG confined in a cigar-shaped trap. These solitons propagate with the velocity 
close to the speed of sound. 



VI. CONCLUSION 



The effective one- dimensional system derived in this paper describes motion of condensate 
in a wide region of parameters when the radial distribution has not reached yet the Thomas- 
Fermi limit and the characteristic axial dimension is not less than the BEG radius in an 
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elongated cigar-shaped trap. This system depends only on one spatial coordinate, what 
can greatly simplify numerical simulations of nonlinear effects in BEC. The advantage of 
our approach is illustrated by applications to description of dark and bright solitons in the 
high density BEC. The one-dimensionality of the effective system permits one to develop 
various analytical approximations. In particular, the dispersion relations for sound and 
radial compression waves are obtained for high density BEC. Evolution of small but finite 
amplitude perturbations is described by the KdV equation which leads to qualitatively new 
behavior of BEC at high enough density: small- amphtude bright solitons can propagate on 
a constant background in a repulsive BEC, whose density is above some threshold value. 
This prediction can be verified by experiment. 

We hope that the quasi-lD theory developed in this paper will find applications in de- 
scription of many other nonlinear processes in the high-density BEC confined in the cigar- 
shaped trap, i.e. when the standard one-dimensional nonlinear Schrodinger equation looses 
its applicability. 
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